APS/123-QED 



Atomic diffraction in counter-propagating Gaussian pulses of laser light 
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We present an analysis of atomic diffraction due to the interaction of an atomic beam with a pair 
of Gaussian light pulses. We derive a simple analytical expression for the populations in different 
diffraction orders. The validity of the obtained solution extends beyond the Raman-Nath regime, 
where the kinetic energy associated with different diffraction peaks is neglected, into the so-called 
channeling regime where accurate analytical expressions have not previously been available for the 
diffraction. Comparison with experimental results and exact numerical solutions demonstrate the 
validity of our analytical formula. 
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I. INTRODUCTION 

Diffraction of light and matter waves lies at the heart 
of a number of diagnostic measurement techniques in 
physics and engineering. Due to their internal level struc- 
ture, atoms and molecules can be diffracted by resonant 
or near resonant laser beams [U |2j [3l HJ |5j [6]. The 
field of atom optics points to applications for high preci- 
sion probing of incrtial effects and fundamental physics 
[3 |9j IHI], as well as to the use of diffraction to di- 
agnose single- and multi-particle properties of the atoms 
themselves in various physical settings [TI] [T3J . The 
use of light induced potentials for structured deposition 
of atoms on surfaces has also been demonstrated [Hj. By 
use of Laguerre-Gaussian laser beams, atomic diffraction 
was recently used to transfer not only linear but also an- 
gular momentum to an atomic gas |151 116] , By tailoring 
the polarization and frequencies of light, it is possible 
to induce complicated multi-level dynamics with various 
potential applications |17j . For most species, using opti- 
cal pumping and suitable polarization schemes, it is also 
possible to restrict the dynamics and to obtain a situa- 
tion of effective propagation of a single component wave 
function in a single periodic potential |17l 118] . 

The coupling of atomic momentum components by ab- 
sorption and stimulated emission events to a comb of 
other momentum states separated by twice the photon 
momentum cannot be solved analytically due to the non- 
commuting kinetic and potential energies. However, for 
two parameter regimes simple analytical solutions have 
been known for a long time. In the Raman-Nath regime 
of short pulse durations, the difference in kinetic en- 
ergy of different momentum states is neglected and the 
coupling leads to an expansion with Bessel-function ex- 
pression for different momentum component amplitudes 
[111 [18]. In the Bragg regime of long and weak pulses, ki- 
netic energy must be strictly conserved by absorption and 
emission of an integer number of counter-propagating 
photon pairs, and the system reduces to an effective, an- 
alytically tractable, two- level system. 

Atomic diffraction is applied extensively in experi- 



ments, and the need for numerical and analytical the- 
oretical approaches beyond the range of validity of the 
Raman-Nath and the Bragg theoretical treatments, is ex- 
plicitly expressed in the recent literature [THJ 120] ■ Here 
we depart from the Raman-Nath approach, and suggest 
an analytical formula which extends into the channel- 
ing regime. Our analytical expression, includes effects 
due to the kinetic energy dispersion, and deploying nu- 
merical simulations we show that it extends the range 
of validity of the simple Raman-Nath approach into the 
channeling regime where both Raman-Nath and Bragg 
descriptions fail. In Section II, we introduce our ana- 
lytical and improved approximate theories. In Sec. Ill 
we present numerical calculations and compare them to 
the different analytical approaches derived in the previ- 
ous section, assessing their validity regime. Finally, Sec. 
IV concludes the paper with a discussion. 



II. THEORY 

We consider an atom with mass M, two atomic levels 
with energy level separation hu>o, and dipole matrix el- 
ement d, interacting with two counter-propagating laser 
field pulses, 

£(z,t) = Si(t- z/c)e^ k "- u ^ +£ 2 {t + z/c)e- i( > kz+UJ V +c.c. 

We imagine that the two pulses arrive at the origin, 
z = 0, symmetrically from opposite directions. For re- 
alistic atomic and field parameters, the width of the 
atomic distribution around z — will be narrow com- 
pared to the spatial range of intensity variation of ei- 
ther light pulse and hence their interaction with the 
atoms is that of a standing wave with a time depen- 
dent strength. We assume that the atoms have initially a 
well defined momentum, po, and we adiabatically elimi- 
nate the excited atomic state, which is off-resonance from 
the laser fields by an amount A = uj — luq, yielding the 
Schrddinger equation for the ground state wave function, 
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FIG. 1: (Color online) Diffraction patterns for various pulse 
durations, r, as predicted by (a) the Raman-Nath solution 
Eq. |BJ|, (b) its modified version, Eq. |6| and (c) the full nu- 
merical solution to the Eq. Q, plotted as functions of the 
momentum, p. The color indicates the population intensity 
such that red corresponds to the maximum and dark blue to 
the minimum. Parameters are: q — 52.9 and po = 20hk. 




FIG. 2: (Color online) Diffraction intensities for a pulse du- 
ration r = 0.0375/(av and intensity q = 52.9 in reduced units. 
The diffraction probabilities are shown as function of the fi- 
nal state momentum. Red diamonds, blue boxes and black 
bullets correspond to the frames (a), (b), and (c) in Fig. [T] 
respectively. 



ip(z, t) = a(z, t) exp(ipoz/h — i-p\tj2MK), 



in—aiz, t) 



where 

VextOM) 



(2) 



2Re(£*(t - z/c)£ 2 (t + z/c)e 



(^(i-z/c^ + l^ + z/c)! 2 



-2ikz 



(3) 



A. Neglecting the variation of the kinetic energy 
for different momentum states 



Neglecting the second order spatial derivative in 
Eq. pj, it becomes solvable by straightforward integra- 
tion |21| . We replace the assumed smooth position de- 
pendence of the amplitude, a(z,t), before the interaction 
with the laser field by a constant, arriving at the solution 

a(z,t) oc exp^l^y" |^| 2 dt' + f |£ 2 | 2 dt' 

J — CO 

+ e -^Hp«t/M-z) J' |£ 1 |2|£ 2 |2 el 2fe P0 t7A/ dt ^^ 

where we will assume real field amplitudes. Physically, 
neglecting the kinetic energy operator amounts to mak- 
ing the dynamics local in position involving only a posi- 
tion dependent phase factor. Assuming Gaussian pulses, 
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£x,2(t±z/c) — £'i ! 2exp{— (t ± z / c) 2 / t 2 } , the definite in- 
tegrals in Eq. Q can be evaluated at t = oo, i.e., after 
the pulse, and the exponential function can be expanded 
in the following series, 



a(z, t) 



i m J 



oo 



a im2k(z-p t/M) 



2v\d\ 2 E l E 2 _ 2{A _y _n 



(5) 



Note that this expression bears resemblance to the 
Raman-Nath solution obtained with constant laser am- 
plitudes [T7J [H]. However, the specific Gaussian shape 
and an exponential dependence on the duration parame- 
ter, t, of the Bessel-function argument highlights certain 
specific physical properties of the diffraction process with 
Gaussian pulses. For all practical purposes, the position 
argument, z, is much smaller than the pulse length cr 
and the first exponential term inside the Bessel-function 
will be replaced by unity in the following. An example 
of the results obtained from Eq. (|5| for the momentum 
distribution, i.e., the absolute value of the square of the 
Bessel-functions evaluated for different diffraction orders, 
m, are plotted in Fig. [ija) as function of the final state 
momentum p = po+2mhk. Note the figure does not show 
the time dependent diffraction, but the asymptotic mo- 
mentum distribution after the pulse, for different values 
of the pulse duration r. 

The Bessel-functions with m =/= all vanish for van- 
ishing arguments, and we see that both short (r — * 0) 
and long (r — * oo) pulses exclude diffraction to new mo- 
mentum components. For short pulses of given strength, 
the explanation is that there is simply no time to signifi- 
cantly diffract the atoms, and the nearly vanishing r ar- 
gument suppresses the higher order Bessel functions. For 
longer pulses the result can be ascribed to the exponential 
factor inside the Bessel-functions, which in turn reflects 
the lack of energy conservation in the diffraction pro- 
cess: In the rest frame of atoms moving at velocity po/M, 
the photons in the counter-propagating light pulses are 
Doppler-shifted in frequency, causing an apparent change 
of 2hkpo/M in the field energy per absorption-emission 
cycle involving both beams. This is a violation of energy 
conservation, permitted because we consider a process of 
finite duration r. Note that the process is suppressed by 
the exponential factor when the corresponding frequency 
mismatch exceeds oc 1/r, in support of our interpreta- 
tion in terms of the energy-time Hciscnberg uncertainty 
relation. Alternatively, the finite duration pulses may be 
thought of as (Gaussian) frequency distributions of light, 
and absorptions and emissions involving photons from 
the counter-propagating pulses may conserve field energy 
in the frame moving at po/M, if the photons involved are 
at frequency components shifted by the amount ±kpo /M 
with respect to the field carrier frequency of the pulses. 
The field power at those frequencies is just given by the 
Gaussian factor, appearing inside the Bessel function ar- 
guments. 



B. Incorporation of the kinetic energy variation 

After having presented the equation governing the 
diffraction within the Raman-Nath approximation, which 
does not take into account the energy difference between 
states with different momenta, it is natural to investi- 
gate if a simple modification of the final result may cap- 
ture the characteristic features of the energy dispersion. 
To this purpose, we shall come back to our interpreta- 
tion of the last exponential factor inside the Bessel func- 
tion in Eq. ([5]). In the previous section we argued that 
this quantity expresses the mean change of field energy 
2hkpo/M per absorption emission cycle, experienced by 
the moving atom due to the Dopper effect, and that it 
is appropriately measured relative to the frequency wit- 
dth of the light pulses due to their finite duration. For a 
Gaussian pulse, the frequency distribution is also Gaus- 
sian, hence the Gaussian dependence inside the Bessel- 
function. Taking into account that the different momen- 
tum states have different kinetic energies, it is thus nat- 
ural to suggest to apply the same argument as above but 
taking into account the change of field energy plus atomic 
kinetic energy. In the frame moving at velocity po/M the 
atom is initially at rest. Scattering process into the mth 
diffraction order by absorption and simulated emission 
of to photons propagating in opposite directions changes 
the field energy by an amount 2mhkp /M. Simultane- 
ously atoms acquire a finite momentum kick of 2mhk in 
the moving frame, and hence their final kinetic energy is 
Am 2 h 2 k 2 /2M. The total energy defect per absorption- 
emission event is therefore 2hk(p + rnhk)/M. Since the 
pulses have Gaussian frequency widths proportional to 
1/r, the atoms may interact with frequency components 
shifted by k(po+mhk)/M in order to conserve total (field 
plus atomic kinetic) energy. Motivated by the above 
analysis we thus suggest to replace po in the exponen- 
tial inside the Bessel-function in Eq. §5§ by po + mhk. 
This modification leads to the expression, 



b( Z) t)= £ 



o im2k(z~p t/M) 



m— — oo 



/2Tr\d\ 2 E 1 E 2 _f(_po±^lJ2z) 2 /2 



h 2 A 



(6) 



for different diffraction orders, m. The probability that 
an atom with initial momentum po acquires the final mo- 
mentum po + 2mhk, is | J m {f{pa + mhk))\ 2 , with 



h 2 A 



(7) 



The results of this modified version of the Raman-Nath 
approximation are shown in Fig. [ljb). We do not claim 
that such alteration is an exact result or even a system- 
atic expansion in any small parameter, and although it is 
physically meaningful it must be justified a posteriori, for 
instance by comparison with an exact calculation. Such 
assessment will be made in the following section. 



III. RESULTS 

In the preceding section we have derived two different 
analytical solutions, which we shall now compare to the 
exact numerical solution of the problem. 



A. Numerical solution 

The full Schrodinger equation without omission of the 
second derivative term is amenable to numerical solution, 
and we have carried out such computations for a wide 
range of atomic and field parameters. Our numerical 
studies enable us to compare the predictions of the an- 
alytical approximations to the exact numerical solutions 
and assess their validity, and in particular, to investigate 
whether the modified form Eq. (JsJ) represent an improve- 
ment over the Eq. d5| or not. To this end, we propagate 
wavef unctions, d(z,t), which interact with the Gaussian 
light pulses, in time from t = — 2r to t = 2r according to 
the Eq. ([2}. The initial wave function represents a parti- 
cle with a mean momentum, p — po, and is modelled as 
a Gaussian wave packet 



d(z, —00) 



_ ~z 2 /2a 2 ipz/h 



(8) 



where the waist, a — 10/fc, is chosen wide enough that 
the momentum spread of the initial wave packet is narrow 
in comparison to the resulting momentum separation of 
adjacent diffraction orders. 

For further convenience, following Rcf. [19 , we intro- 
duce dimcnsionless parameters for time and laser field 
intensity 



d 2 E l E 2 



Ah 2 uj r A 



u> r t, 



(9) 



where u> r — hk 2 /2M is the photon recoil frequency. 



B. Comparison of different analytical solutions 

In Fig[lja)-(c), we have plotted the diffraction patterns 
for a range of pulse durations, r, as functions of the mo- 
mentum, p for two different approximations and the ex- 
act solution to the diffraction problem, respectively. We 
obtain the diffraction peak intensities as a square mod- 
ulus of the Fourier transformed wave functions. Note 
that the 'onion shape' does not show the momentum dis- 
tribution as a function of time, instead, each horizontal 
cut through the figure indicates the asymptotic momen- 
tum distribution after a pulse of the given duration. The 
pulse parameters for all of the frames are q = 52.9 and 
Po = 20Hk. The modified Raman-Nath solution, shown 
in Fig. [ljb) , is observed to be in an excellent agreement 
with the exact numerical solution, Fig. [ijc), for these 
parameters . 

To obtain better quantitative feeling for the quality of 
the equation ([6]), we plot the diffraction pattern for a 




FIG. 3: (Color online) Diffraction patterns for various pulse 
durations, t, as predicted by the modified Raman-Nath solu- 
tion Eq. Q, plotted as a function of the momentum, p. The 
color indicates the amplitude such that red corresponds to the 
maximum and dark blue to the minimum. The pulse intensity 
q = 3.5 and the atom arrives to the interaction region with 
an initial momentum p = hk. 




FIG. 4: (Color online) Diffraction patterns for various pulse 
durations, r, as predicted by the full numerical solution to 
Eq. j2j, plotted as a function of the momentum, p. The color 
indicates the amplitude such that red corresponds to the max- 
imum and dark blue to the minimum. The pulse intensity 
q = 3.5 and the atom arrives to the interaction region with 
an initial momentum p — hk. 



single pulse duration, r = 0.0375/uv as shown in Fig. [2] 
Here the red diamonds, blue boxes and black bullets cor- 
respond to the values shown in the frames (a), (b), and 
(c) in Fig. [l] respectively. From this graph, it is clear 
that our modified solution still works excellently in the 
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pulse duration [r T 



marizes the applicability of various approximations on 
the (r r , q) parameter space. In the spirit of Ref. |19J . the 
region boundaries are defined by r r — \/2^/2q and q = 1. 
Black dots mark the experiments as in Ref. [19 . The 
upper horizontal line connecting triangle symbols corre- 
sponds to the results shown in Fig. [ijb) and the lower line 
connecting square symbols depicts the parameters used 
in Fig. [3] The dashed part of the lower line indicates that 
Eq. (|6| fails for these parameters, deep into the channel- 
ing regime. However, it is clear that the validity of our 
analytical solution enters the previously unknown chan- 
neling territory. It is also worth noting that for the case 
of high intensity (high q value), demonstrated in Fig. [T| 
our expression maintains its validity for all pulse dura- 
tions, r, since the diffraction process itself is confined to 
relatively short times. 



FIG. 5: (Color online) Pulse duration versus intensity dia- 
gram for atomic diffraction dividing the parameter space in 
three distinctively different regimes. Black bullets correspond 
to the experiments in Ref. [19] . Triangles joined by a line cor- 
respond to the modified Raman-Nath solution in Fig. [ljb) 
and the boxes joined by a line are from Fig. [3] Dashed line 
marks the region where the modified Raman-Nath formula 
loses its validity. 



regime where the Raman-Nath approximation fails. 



C. Validity assessment and relation to experiments 

We have also compared our modified Raman-Nath pre- 
diction, Eq. Q, with the actual experiment as reported 
in Ref. [T!5] in addition to the full numerical solution. 
In the experiment, the smooth-shaped pulse was varied 
as cos 2 (at) which is well approximated by our Gaussian 
shaped pulse and allows a direct comparison. With pa- 
rameters corresponding to the experiment, q — 3.5, r r = 
[0,1.5] and po = hk, Fig. [3] shows our prediction while 
Fig. [4] is the corresponding numerical solution. Compari- 
son between these and the experimental result presented 
in jT9] , we observe all of them to be in fair agreement for 
short pulse durations, i.e., in the Raman-Nath regime. 
For longer pulse durations, both experiment and numeri- 
cal solution exhibits Bragg-like channeling behavior with 
an associated Rabi cycling. Such long-time behavior for 
these parameters cannot be reproduced by our analytical 
formulas. This may be because the experiment is car- 
ried out with parameters very far from the Raman-Nath 
regime. The experiment is in fact not that far from the 
Bragg regime, where we know that a correct treatment 
of the kinetic energy is crucial. 

Figure p] is in parts reproduced from Ref. [19] and sum- 



IV. DISCUSSION 

We have presented an analytical solution to a problem 
of atomic diffraction from a pulsed time-dependent light 
grating for the so called Raman-Nath regime where the 
kinetic energy of the atoms participating in the diffrac- 
tion process may be neglected. We have extended this 
approximation by including certain aspects of the kinetic 
energy in the diffraction process and we have proposed 
a phcnomcnological improvement to this solution which 
is based on physically motivated arguments. Through 
a comparison with an exact numerical solution of the 
Schrodinger equation, it is shown that our modification 
extends the validity regime of the parameter space of 
the analytical solution from the Raman-Nath into the 
channeling regime. Qualitatively the modification pre- 
dicts the asymmetry due to the Doppler effect in the 
diffraction pattern. The proposed analytical formula is 
compared with available numerical results and is found 
to be in good agreement even fairly deep in the chan- 
neling regime of strong and long interactions. Although, 
containing certain elements of Bragg diffraction, the pro- 
posed equation is not, however, capable of predicting 
correctly the effective two-state Rabi cycling between 
momentum states ±p . Exploration of the channeling 
regime from the Bragg-formalism seems a promising av- 
enue for attempts to provide good analytical expressions 
for atomic diffraction closer to the Bragg regime. To- 
gether with the present work, this may perhaps exhaust 
the full parameter range with good analytical approxi- 
mations for the atomic diffraction. 

The authors acknowledge support from a NATO link- 
age grant, and A. M. acknowledges support via Yerevan 
State University grant NFSAT/CRDF UCEP-02/07. 



[1] Phillip L. Gould, George A. Ruff, and David E. [2] O. Carnal and J. Mlynek, Phys. Rev. Lett. 66, 2689 
Pritchard, Phys. Rev. Lett. 56, 827 - 830 (1986). 



(1991). 

[3] D.W. Keith, C. R. Ekstrom, Q. A. Turchette, and D. E. 

Pritchard, Phys. Rev. Lett. 66, 2693 (1991). 
[4] F. Riehle, Th. Kister, A. Witte, J. Helmcke, and Ch. 

Borde, Phys. Rev. Lett. 67, 177 (1991). 
[5] F. Shimizu, K. Shimizu, and H. Takuma, Phys. Rev. A 

46, R17 (1992). 
[6] M. Kasevich and S. Chu, Phys. Rev. Lett. 67, 181 (1991). 
[7] J. M. McGuirk, G. T. Foster, J.B. Fixler, M. J. Snadden, 

and M. A. Kasevich, Phys. Rev. A 65, 033608 (2002). 
[8] T. L. Gustavson, P. Bouyer, and M. A. Kasevich, Phys. 

Rev. Lett. 78, 2046 (1997). 
[9] A. Lenef, T. D. Hammond, E. T. Smith, M. S. Chapman, 

R. A. Rubenstein, and D. E. Pritchard, Phys. Rev. Lett. 

78, 760 (1997). 

[10] B. Canuel, F. Leduc, D. Holleville, A. Gauguet, J. Fils, A. 

Virdis, A. Clairon, N. Dimarcq, Ch. J. Borde, A. Landra- 

gin, and P. Bouyer, Phys. Rev. Lett. 97, 010402 (2006). 
[11] J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper- 

Kurn, D. E. Pritchard, and W. Ketterle, Phys. Rev. Lett. 

82, 4569 (1999). 
[12] Yu. B. Ovchinnikov, J. H. Miiller, M. R. Doery, E. J. 

D. Vredenbregt, K. Helmerson, S. L. Rolston, and W. D. 

Phillips, Phys. Rev. Lett. 83, 284 (1999). 
[13] J. Steinhauer, R. Ozeri, N. Katz, and N. Davidson, Phys. 

Rev. Lett. 88, 120407 (2002). 
[14] B. Rohwedder, Am. J. Phys. 75, 394 (2007). 
[15] M. F. Andersen, C. Ryu, P. Clade, V. Natarajan, A. 

Vaziri, K. Helmerson and W. D. Phillips, Phys. Rev. Lett. 

97, 170406 (2006). 
[16] T. P. Simula, N. Nygaard, S. X. Hu, L. A. Collins, B. I. 

Schneider, and K. M0lmer, cond-mat/ '1234.5678 (2007). 
[17] C. S. Adams, M. Siegel, and J. Mlynek. Physics Reports, 

240, 143 (1994). 
[18] Atom Interferometry, edited by R. Paul Berman (Aca- 
demic, London, 1997). 
[19] C. Keller, J. Schmiedmayer, A. Zeilinger, T. Nonn, S. 

Diirr, and G. Rempe, Appl.Phys. B 69, 303 (1999). 
[20] H. Miiller, S. Chiow, and S. Chu, |arXiv:0704.2627| j(2007). 
[21] V.M.Arutunyan and A.Zh.Muradyan, Dokladi Akad. 

Nauk Arm. SSR, 60, n.5, 275 (1975). 



